Mon. Not. R. Astron. Soc. 000. HH291 (2010') Printed 29 September 2010 (MN M£X style file v2.2) 



Multi-wavelength polarimetric study towards the open 
cluster NGC 1893 

o : 
o 

C. Eswaraiah, 1 * A. K. Pandey, 1 G. Maheswar, 1 Biman J. Medhi, 1 
OO ; J. C. Pandey 1 , D. K. Ojha 2 and W. P. Chen 3 



00 



ON 

o 
o 



1 Aryabhatta Research Institute of Observational Sciences, Manora Peak, Nainital 263129, India 

2 Tata Institute of Fundamental Research, Mumbai 400 005, India 
3 Institute of Astronomy, National Central University, Chung-Li 32054, Taiwan 

00 ' 

Oh' Accepted , Received ; in original form 

6 

ABSTRACT 

. We present multi-wavelength linear polarimetric observations for 44 stars of 

OO 



the NGC 1893 young open cluster region along with y-band polarimetric 



' observations of stars of other four open clusters located between I ~ 160° 

in : 

to 175°. We found evidence for the presence of two dust layers located at a 



distance of ~ 170 pc and ~ 360 pc. The dust layers produce a polarization 
Py ~ 2.2%. It is evident from the clusters studied in the present work that, in 
the Galactic longitude range I ~ 160° to 175° and within the Galactic plane 
(|fe|<2°), the polarization angles remain almost constant, with a mean ~ 163° 
and a dispersion of 6°. The small dispersion in polarization angle could be 
due to the presence of uniform dust layer beyond 1 kpc. Present observations 
reveal that in case of NGC 1893, the foreground two dust layers, in addition to 
the intracluster medium, seems to be responsible for the polarization effects. 
It is also found that towards the direction of NGC 1893, the dust layer that 
exists between 2-3 kpc has a negligible contribution towards the total observed 
polarization. The weighted mean for percentage of polarization (P max ) and 
the wavelength at maximum polarization (\ m ax) are found to be 2.59 ±0.02% 
and 0.55 ± 0.01 /im respectively. The estimated mean value of X m ax indicates 
that the average size of the dust grains within the cluster is similar to that 
in the general interstellar medium. The spatial variation of the polarization is 
found to decrease towards the outer region of the cluster. In the present work, 
we support the notion, as already has been shown in previous studies, that 
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polarimctry, in combination with (U — B) / (B — V) colour-colour diagram, is 
a useful tool for identifying non-members in a cluster. 

Key words: polarization - dust, extinction - open clusters and associations: 
individual: (NGC 2281, NGC 1664, NGC 1960, Stock 8 and NGC 1893). 



1 INTRODUCTION 

Dust properties govern several physical and chemical phenomena in the interstellar medium 
(ISM) and act as a tracer of local environmental conditions. When starlight passes through 
various components of interstellar dust grains, which are aligned to the Galactic magnetic 
field, the radiation becomes partially plane-polarized, typically at the level of a few percent. 
The nature of this polarization reveals important information regarding the shape, size 
and composition of interstellar dust. The detailed process(es) by which grains get aligned 
with the magnetic fields has (have) long been, and is still actively under study (Davis & 
Greenstein 1951; Dolginov 1990; Cho & Lazarian 2005). In general the dust grains tend to 
align with their long axis perpendicular to the magnetic field (Purcell 1979). The clue to 
the interstellar origin of the polarization came from the observed correlation between the 
degree of polarization (P) and the colour excess, E(B — V). The values of P (in visual 
wavelengths) for the stars with large E(B — V) are found to be in the range between zero 
and a maximum value given by P max /E(E> — V)= 9% mag^ 1 (Aannestad & Purcell 1973). 
The relation between P ma x and the colour excess, and the variation of P with wavelength, 
are interpreted in terms of the grain properties and the efficiency of the grain alignment. 
Therefore, polarimetry is a useful technique to investigate the properties like maximum 
polarization P\ max , the wavelength X max corresponding to P\ max and the orientation of the 
magnetic field in various Galactic locations. 

The photometric and spectroscopic information already available for stars in open clus- 
ters are specifically important to make a meaningful study of the dust grains located in the 
foreground and in the intracluster regions (e.g., Trumpler 27, Feinstein et al. 2000; NGC 
6231, Feinstein et al. 2003; Stock 16, Feinstein et al. 2003; Hogg 22 & NGC 6204, Martinez 
et al. 2004; NGC 5606, Orsatti et al. 2007; NGC 5749, Vergne et al. 2007; IC 1805, Medhi et 
al. 2007; NGC 6250, Feinstein et al. 2008; NGC 654, Medhi et al. 2008; NGC 6124, Vergne 
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Table 1. The basic parameters of the clusters. 



Cluster Id 


1 


b 


Distance 


E(B-V) 


log(age) 


Reference 




(°) 


(°) 


(pc) 


{mag) 


(yr) 




NGC 2281 


174.90 


+16.88 


558 


0.06, 0.11 


8.70 


Karchenko et al. (2005), Glaspey (1987) 


NGC 1664 


161.68 


-0.45 


1199 


0.25 


8.72 


Karchenko et al. (2005) 


NGC 1960 


174.54 


+1.07 


1330 


0.22 


7.40 


Sharma et al. (2006) 


Stock 8 


173.37 


-0.18 


2050 


0.40-0.60 


6.00-6.70 


Jose et al. (2008) 


NGC 1893 


173.59 


-1.68 


3250 


0.40-0.60 


6.60 


Sharma et al. (2007) 



Table 2. The observational details 



Cluster Id 


Date of Observation 
(year, month, date) 


Passband(s) 


NGC 2281 


2009, 12, 24 


V 


NGC 1664 


2009, 12, 24 


V 


NGC 1960 


2009, 11, 23 


V 


Stock 8 


2009, 11, 23 


V 


NGC 1893 


2008, 11, 8 & 9 


BV(RI) C 



et al. 2010; NGC 6823, Medhi et al. 2010). In this paper, we present multi- wavelength po- 
larimetry of stars towards NGC 1893. We also present l/-band polarimetry of stars towards 
four additional clusters, viz NGC 2281, NGC 1664, NGC 1960 and Stock 8, to get informa- 
tion about the ISM foreground to the cluster NGC 1893. The main aim of this study is to 
investigate the dust properties as a function of distance towards the anticenter direction of 
the Galaxy (I ~ 160° to 175°) using stars in open clusters. 

The basic parameters of the observed clusters are given in Table [TJ Using these five 
open clusters, we have made an attempt to study the dust properties of the ISM distributed 
between ~ 600 pc (NGC 2281) and ~ 3 kpc (NGC 1893). The paper is organized in the 
following manner. In section |2j we present a brief discussion on the observations and data 
reduction. We present our results in section |3l and analysis and discussion in section HI The 
dust components responsible for the observed polarization are discussed in section [5] We 
conclude our results in section |6l 

2 OBSERVATIONS AND DATA REDUCTION 

Polarimetric observations were carried out using the ARIES Imaging Polarimeter (AIMPOL: 
Rautela, Joshi & Pandey, 2004; Medhi et al. 2007, 2010) mounted at the Cassegrain focus of 
the 104-cm Sampurnanand telescope of the Aryabhatta Research Institute of observational 
sciencES (ARIES), Nainital, India, coupled with a TK 1024x1024 pixel 2 CCD camera. 
The AIMPOL consists of a half-wave plate (HWP) modulator and a Wollaston prism beam- 
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splitter. The observations were carried out in B, V, R c and I c (As ..=0.440 fim, Ay e// =0.530 
jUin, Ai? Ce// =0.670 /im and Aj ..=0.800 /im) photometric bands. Details of the observations 
are given in Table [2 Each pixel of the CCD corresponds to 1.73 arcsec and the field-of-view 
is ~ 8 arcmin in diameter. The FWHM of the stellar image varied from 2 to 3 pixels. The 
read out noise and gain of the CCD are 7.0 e _1 and 11.98 e _1 /ADU, respectively. Due to 
the absence of a grid in AIMPOL, we manually checked for any overlap of ordinary and 
extra-ordinary images of the sources. 

Fluxes of ordinary (J ) and extra-ordinary (I e ) beams for all the observed sources with 
good signal-to-noise ratio were extracted by standard aperture photometry after bias sub- 
traction using the iraJ^ package. The ratio R(a) is given by: 
W _ 1 

R ( a ) = f|j— 7 = Pcos ( 2 ^ - 4 «)> C 1 ) 

where P is the fraction of the total linearly polarized light and 9 is the polarization angle of 
the plane of polarization. Here a is the position of the fast axis of the HWP at 0°, 22.5°, 45° 
and 67.5° corresponding to the four normalized Stokes parameters respectively, q [R(0°)], 
u [R(22.5°)], qi [R(45°)] and U\ [R(67.5°)]. The detailed procedures used to estimate the 
polarization and polarization angles for the programme stars are described by Ramaprakash 
et al. (1998), Rautela, Pandey, & Joshi (2004) and Medhi et al. (2010). Since polarization 
accuracy is, in principle, limited by photon statistics, we estimated the errors in normalized 
Stokes parameters <jR( a ) (<J q , a u , a qi and a ui in %) using the expression (Ramaprakash et 
al. 1998): 



VR(«) = V(^ e + N + 2N b )/(N e + N ) (2) 

where N e and N are the counts in extra-ordinary and ordinary rays respectively, and Nj,[= 
N be +N bo -^ - g ^ e avera g e background counts around the extra-ordinary and ordinary rays of 
a source. The individual errors associated with the four values of R(a), estimated using 
equation (j2J), are used as weights in calculation of P and 9 for the programme stars. 

To correct the measurements for null polarization (or instrumental polarization) and the 
zero-point polarization angle, we observed several polarized and unpolarized standards taken 
from Schmidt, Elston & Lupie (1992). The results on polarized and unpolarized standards 
are given in Table [3J The values of 9 are in equatorial coordinate system measured from 
the North increasing towards the East. Both the observed degree of polarization [P(%)] 

1 IRAF is distributed by National Optical Astronomical Observatories, USA. 
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Figure 1. Stellar polarization superimposed on the 20 X 20 i2-band DSS II image of the field containing NGC 1893. The 
length of each polarization vector is proportional to Py. A vector with a polarization of 2% is drawn for reference. The dash- 
dotted line marks the Galactic parallel 6=-1.53°. Ha emission stars identified by Marco and Negueruela (2002) (S2R2N43) and 
by Negueruela et al (2007) (S1R2N35, S1R2N56 and S1R2N44) are shown with cross symbols. The identification numbers are 
taken from Cuffey and Shapley (1937). The circle represents the boundary (~ 6 radius, Sharma et al. 2007) of the cluster. 



and polarization angle [0(°)} for the polarized standards are in good agreement with those 
given by Schmidt et al. 1992. The observed normalized Stokes parameters q and u (q%, u%) 
for standard unpolarized stars are also given in Table [3j The instrumental polarization of 
AIMPOL on the 104-cm Sampurnanand Telescope has been monitored since 2004 in different 
projects and found to be less than 0.1% in different bands (Rautela, Pandey, & Joshi, 2004; 
Medhi et al. 2007, 2008, 2010; Pandey et al. 2009). 
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3 RESULTS 

3.1 NGC 1893 

The polarization measurements for stars in the region of NGC 1893 are listed in Table |H 
The star identification numbers (column 1) are taken from Cuffey and Shapley (1937). The 
V^-band magnitudes are given in column 2. The degree of polarization P (in percent) and 
polarization angles 9 (in degree) measured in B, V, (R, I) c bands and their corresponding 
standard errors (ep and eg) are given in columns 3 to 10. 

The sky projection of the l/-band polarization vectors for the 44 stars measured towards 
the NGC 1893 region is drawn on the i?-band Digitized Sky Survey II (DSS II) image (Fig. 
[Q). The length of each polarization vector is proportional to the degree of polarization. A 
vector with a polarization of 2% is drawn for reference. The dash-dotted line superimposed 
in Fig. [T] is indicating the orientation of the projection of the Galactic plane (GP) which has 
a polarization angle of Qqp = 141°. The average polarization (P) and polarization angle (9) 
for 44 stars are found to be 2.4 ± 0.6%, 160 ± 5° in B, 2.6 ± 0.7%, 162 ± 5° in V, 2.5 ± 0.7%, 
159 ± 6° in R c and 2.2 ± 0.5%, 161 ± 6° in J c bands, respectively. 

3.2 NGC 2281, NGC 1664, NGC 1960 and Stock 8 

Polarimetric measurements in ^/-band for 14 stars towards NGC 2281, 27 stars towards NGC 
1664, 15 stars towards NGC 1960 and 21 stars towards Stock 8, have also been carried out. 
The results are listed in Tables [51 [7] and |8] respectively. The star identification numbers for 
NGC 2281, NGC 1664, NGC 1960 and Stock 8 are taken from Vasileviskis & Balz (1959), 
Larsson-Leander (1957), Boden (1951) and Mayer (1964) respectively. 

The mean values of P v and 9 V towards NGC 2281, NGC 1664, NGC 1960 and Stock 8 
are found to be 0.9 ± 0.2% and 16 ± 6°; 1.8 ± 0.7% and 172 ± 8°; 1.2 ± 0.2% and 159 ± 5° 
and 2.4 ± 0.6% & 160 ± 7°, respectively. 

4 ANALYSIS AND DISCUSSION 

Fig. [2] presents the distribution of Py and 6y for stars in the observed cluster regions. Four 
clusters studied here, namely NGC 1664, NGC 1960, Stock 8 and NGC 1893 are located at 
Galactic latitude |6| < 2°. The nearest open cluster NGC 2281 studied here is located at 
b = +16.88°. A Gaussian fit to the distribution of Py and Qy yields a mean and a standard 
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Figure 2. Histograms for P v (left panel) and 6 V (right panel) in five observed clusters, NGC 2281, NGC 1664, NGC 1960, 
Stock 8 and NGC 1893. The polarization angle for the Galactic plane of each cluster is shown with a continuous line and is 
labeled with GP. Each histogram is fitted with a Gaussian function. 



deviation as 0.8 ± 0.2% and 15 ± 9° for NGC 2281; 1.3 ± 0.2% and 170 ± 5° for NGC 1664; 
1.2 ± 0.2% and 159 ± 3° for NGC 1960; 2.3 ± 0.1% and 160 ± 3° for Stock 8; and 2.6 ± 0.6% 
and 160 ± 4° for NGC 1893. It is evident from Fig. [2] that the mean value of Py increase 
with the distance of a cluster. 

The polarization angle corresponding to the Galactic plane for each cluster is drawn with 
a continuous line. Fig. [2] (right panel) indicates that the difference between GP and mean 
9 V is the largest for the nearby least polarized (0.8 ± 0.2%) cluster (NGC 2281: Z=174.90° 
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and 6=+16.88°). Whereas, in the case of four clusters, namely NGC 1664, NGC 1960, Stock 
8 and NGC 1893, the difference between GP and mean Qy (Gaussian) is small and is almost 
the same. Fig. [2] also reveals that in the Galactic longitude range I ~ 160° to ~ 175° and 
within the Galactic plane (|6| < 2°), the polarization angles remain almost constant. The 
average of the mean values of polarization angles for the four clusters is found to be 163° 
with a standard deviation of 6°. This implies that at least beyond ~1 kpc the magnetic 
field orientation remains almost unchanged. However, for a smaller distance (~ 0.6 kpc), 
the polarization angle towards the northern hemisphere (NGC 2281) of the Galactic plane 
is found to be 16° with a standard deviation of 6°. 

The polarization towards the Taurus molecular complex (I = 174.13°, b = —13.45°), 
taken from Heiles (2000), reveals a mean polarization angle of 58° with a standard deviation 
of 38°. The region selected around the Taurus molecular complex is bound between the 
Galactic longitude I = 160° to 175° and latitude b < —2°. It is interesting to mention that, 
the mean polarization angles at high Galactic latitude > 2°) are significantly different 
from those in the Galactic plane |6| < 2°. 

4.1 Member identification 

The individual Stokes parameters of the polarization vector of the \^-band, Py, given by 
Qv = Py cos(26V) and Uy = Py sm(29 v ) are estimated for all the observed stars towards 
five clusters and presented on a Uy versus Qy plot in Fig. [3j In such a diagram, members 
of a star cluster are expected to group together, while non-members are expected to show 
a scattered distribution. The reason is that the measured polarization of a star depends on 
the column density of aligned dust grains that lie in front of the star and hence the degree 
of polarization would be similar, lower or higher depending on whether the star is a cluster 
member, a foreground or a background to the cluster. Likewise, the polarization angles of 
cluster members would be similar but could be different for foreground or background non- 
member stars as light from them could have contributions from different or additional dust 
components. Therefore Uy — Qy plot could be a useful tool to segregate the members and 
non-members of a cluster. Stars with intrinsic polarization, e.g., due to an asymmetric distri- 
bution of matter around young stellar objects (YSOs) and/or rotation in their polarization 
angles (see section 14. 2[) may also create scattered distribution in the Uy — Qy plane. 

As shown in Fig. [3j in the case of nearby clusters viz NGC 2281, NGC 1664 and NGC 
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Figure 3. The U v versus Q v diagrams for stars towards NGC 2281, NGC 1664, NGC 1960, Stock 8 and NGC 1893 cluster 
regions, measured in the V-band. In each case a box with dashed line marks the boundary of mean Py ± a and mean 0y ± a. 
Stars lying outside the la box are shown by filled circles. It can be seen that the majority of the stars of apparent grouping lie 
within the box (open circles). The apparent grouping in the case of NGC 2281, NGC 1664 and NGC 1960 is bound by a dotted 
circle. The thin continuous curve in the case of Stock 8 and NGC 1893 separates foreground field stars from cluster members. 
The dotted curve in the case of NGC 1893 demarcates the stars having relatively high polarization and high E(B — V) values 
as well as having intrinsic polarization (see the text). 
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Figure 4. (U — B) versus (B — V) colour-colour diagram for the stars in the cluster region. Theoretical zero-age-main sequence 
(ZAMS) is taken from Schmidt-Kaler (1982) and is shifted along a reddening vector with an adopted slope of E(U — B) /E(B — 
V)=0.72, to match the observed colours. The symbols are same as in Fig. E] In the case of NGC 1893, the Ha emission stars are 
indicated with cross symbols. The lowermost-right panel shows the enlarged view of the colour-colour diagram shown within 
the inset of NGC 1893 colour-colour diagram (lower left panel). 



1960, a grouping (bound visually by the dotted circle) is apparent, whereas in the case 
of distant clusters i.e. Stock 8 and NGC 1893, the Uy — Qy diagram shows a scattered 
distribution. A scattered distribution in the case of a distant cluster is expected as the 
cluster region is contaminated by foreground/ background field stars. Moreover, clusters 
Stock 8 and NGC 1893 are young and have differential reddening, hence cluster matter may 
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also affect the polarization. The apparent clustering in the case of NGC 2281, NGC 1664 
and NGC 1960 further supports the above notion as these are relatively nearby clusters and 
have less field star contamination. Since these clusters are of intermediate age (> 25 Myr), 
the effect of variable reddening is absent. The distribution is significantly scattered in the 
case of NGC 1893, which is the most distant cluster in the present study. It seems difficult 
to identify probable cluster members of NGC 1893 using the Uy — Qy diagram as shown in 
Fig. |3) To further elucidate the membership in the clusters which do not show any apparent 
grouping we plot a box with dashed line in Uy — Qy plots (Fig. [3]) having boundaries of mean 
Py ±0" and mean 8y±a (as mentioned in section E])- It can be seen that the majority of stars 
of apparent grouping lie within the boundaries of the box as mentioned above. Therefore, 
in the case of clusters, where grouping is absent (Stock 8 and NGC 1893), we presume that 
the stars lying within the lcr box of mean Py and Qy values are probable members of the 
cluster. 

The combination of Uy — Qy diagram and (U — B)/(B — V) colour-colour diagram can 
yield a better identification of probable members in a distant cluster (see e.g. Haikala 1995, 
Feinstein et al. 2008). Fig. H] shows the (U — B)/(B — V) colour-colour diagram for NGC 
2281, NGC 1664, NGC 1960, Stock 8 and NGC 1893 cluster regions. For the clusters NGC 
2281, NGC 1664, NGC 1960 and Stock 8, the (B - V) and (U - B) colours are taken from 
the Tables El El [7] and El respectively, whereas for the cluster NGC 1893, the colours are 
taken from Sharma et al. (2007) and Massey et al. (1995). The continuous curve represents 
the reddened ZAMS for the cluster region as per E(B — V) values mentioned in Tabled! In 
Fig. El the stars lying within the la box are shown by open circles, whereas stars outside the 
la box are shown by filled circles. It is apparent from Fig. El that stars located within the 
la box follow the general reddening of the cluster region, hence may be probable members 
of the clusters. On the basis of Figs El and HI we can draw the following conclusions: 

NGC 2281: All the stars lying within the la box and the boundary of apparent grouping 
follow the ZAMS reddened for E(B — V) (=0.11 mag) of the cluster region. Although stars 
#55, 56 and 74 lie on the reddened ZAMS, they have (B — V) > 1.0. These stars could not 
be MS members of the cluster having log(age) = 8.7 yr. The star #60 is located away from 
the grouping (Fig. E]) • However it follows the reddened ZAMS and its colours are comparable 
to the probable MS members. 

NGC 1664: Stars located within the region bound by a dotted circle (see Fig. E]) have Py 
< 2% and nicely follow the ZAMS reddened by E(B — V) = 0.25 mag. Hence these could be 
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probable members of the cluster. The stars #1, 22, 27 and 79 are located outside the circle 
but within the la box. These stars show polarization in the range 2.0 < Py < 2.5% and 
follow the ZAMS reddened by E(B - V) ~ 0.55 mag. Stars #9, 13, 47 and 87 lie outside the 
la box. These stars show Py values ranging from 2.5 — 4.1 %. These two groups of stars 
should be background stars. Star #13 has Py = 2.59%, but its estimated E(B — V) ~ 0.25 
is comparable to the cluster's E(B — V) value. As the age of the cluster is ~ 500 Myr, the 
star #13 (V— 14.36, 5- V=1.09) could not be a MS member. Star #164, located significantly 
away from the grouping, is the lowest polarized star (0.97%) and its polarization angle is 
significantly different (26±7°) from those of the remaining stars of the region. It should be 
a foreground non-member. 

NGC 1960: The stars located within the bound region and the la box indicate E(B— V) ~ 
0.22 mag. Star #55 is located significantly away from the general distribution, but follows the 
ZAMS reddened by the average reddening of the region (cf. Fig. H]). This star has relatively 
a high polarization (Py = 1.44 ± 0.24%) and small polarization angle (9y = 147 ± 5°) 
in comparison to the remaining stars of the region. It should be a background field star. 
Star #61, the least polarized star (Py = 0.81%) of the region, lies on the reddened ZAMS, 
but is located away from the general distribution. We conclude that star #61 is probably 
a foreground field star. Although star #77 is located beyond the la box, it follows the 
reddened ZAMS. This star has polarization Py = 1.50 ± 0.31% and could be a background 
non-member. Star #91 seems to be associated with the apparent grouping. Its polarization 
(Py = 1.25 ±0.14%) is comparable to the cluster mean Py value, hence it could be a member 
of the cluster. 

Stock 8: The Uy — Qy diagram does not show any grouping. The (U — B)/(B — V) 
colour-colour diagram indicates that barring stars #211 and 365, the remaining stars lying 
within the la box have E(B — V) in the range of 0.40 — 0.60 mag, hence should be members 
of the cluster. Stars #11, 160 and 173 lie outside the la box. These three stars, along with 
stars #19, 211 and 365, have E(B — V) ~ 0.25 mag. Hence these should be foreground 
stars. Stars #16 and 231 are distributed near the edge of the la box and have E(B — V) 
and polarization values comparable to the mean values for the cluster region. These two 
stars could be cluster members. The stars #275 and 325 are located out side the la box, 
however they seem to be reddened B type stars of the region (cf. Fig. H|). Star #275 has 
significantly different polarization angle (136±5°) from the mean value for the cluster region 
(160±7°), although the polarization value (2.37±0.39%) is similar to the cluster mean value 
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(2.40±0.60%). Star #325 has a relatively higher P v (4.18±0.44%) value than the mean value 
for the cluster region (2.40±0.60%) although the polarization angle (165±3°) is similar to 
the cluster mean value (160±7°). Hence the membership determination of these two stars, 
#275 and 325, is uncertain. The approximate boundary to demarcate the foreground stars 
from the cluster members is shown by a thin continuous curve (see Fig. [3]). 

NGC 1893: The colour-colour diagram of the region reveals that barring stars #197 and 
356, all the stars lying within the la box have EiB — V) ~ 0.4—0.6 mag, indicating that 
these stars may be cluster members. Fig. @] reveals that stars #48, 72, 73, 92, 209, 310, 
327, 329 and 355 have E(B — V) ~ 0.2 mag, hence these should be field stars. Five stars 
#33, 134, 168, 219 and 222 are located near the boundary of the la box. These stars have 
E(B — V) in the range of 0.40—0.60 mag, hence these could also be members of the cluster. 
Star #149 is apparently located away from the distribution. This star has relatively higher 
value of e (2.43; cf. Sec. 4.2 and Table [9]) and may have a rotation in the polarization angle. 
Two stars #141 (S2R2N43, Marco and Negueruela 2002) and #196 (S1R2N35, Negueruela 
et al. 2007) are Ha emission stars. On the basis of polarization, reddening and position in 
the Uy — Qy plot, star #141 should be a cluster member. However, star #196 is a highly 
reddened MS star (E(B — V)=0.83) and is probably an HAe/Be star (see Fig. E]), but its 
Py (2.45 ± 0.29%) and 9 V (160 ± 3°) suggest that it could have undergone depolarization 
effect. Hence, the membership of star #196 is uncertain. The stars having relatively high 
polarization and high E(B — V) values as well as having intrinsic component of polarization 
have been demarcated using a dotted curve as shown in Fig. [3j Fig. H] indicates that stars 
#82 and 98 are highly reddened stars, whereas 123 and 199 are Ha emission stars with 
possible intrinsic component of polarization. Star #108 seems to be a classical Be star (cf. 
Fig. E]) and may have an intrinsic component of polarization. Star #213 may have rotation 
in its polarization angle as its e = 1.45 (see Table [9]). 

The probable members of the cluster identified using Uy — Qy and colour-colour dia- 
grams, along with the kinematic membership probability (Mp) taken from Vasileviskis & 
Balz (1959) for NGC 2281 and from Dias et al. (2006) for remaining clusters, are mentioned 
in Tables El El [TJ M and [9j Stars with kinematic membership probability ^ 50% are con- 
sidered as members. The member and non-member stars are represented with M and NM 
respectively. Stars with uncertainty in their membership determination are indicated with a 
'?' symbol. A comparison indicates that the discrepancy between the membership estimated 
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from kinematic criterion and in the present work is ~ 15% — 20% for the nearby clusters (i.e., 
NGC 2281, NGC 1664 and NGC 1960). The discrepancy is found to increase (~ 30% -45%) 
for distant clusters (i.e., Stock 8 and NGC 1893). The mean values of Py and Qy of member 
stars of the respective clusters as identified above as well as of field stars lying in the cluster 
regions are given in Table [TOj The mean values of Py for members of the two nearby clusters 
viz NGC 2281 and NGC 1960 are the same as those for field stars in these regions, whereas 
in the case of other nearby cluster (NGC 1664), the non- members show higher value for 
mean Py. The higher values of Py for field stars towards the direction of cluster NGC 1664 
are due to the fact that the non-member field stars are located in the background of the 
cluster and have higher extinction (see Fig. H|). In the case of the distant clusters Stock 8 
and NGC 1893, the mean values of Py for non-members are lower than those for member 
stars. The mean values of Py obtained for cluster members further manifest an increasing 
trend with distance. Table ITOl also indicates that the mean 9y values are almost the same for 
members and non-members lying towards the direction of the clusters studied in the present 
work. 

The Uy — Qy plot is also a useful tool to study the interstellar dust distribution as a 
function of distance from the Sun to the clusters (e.g., Feinstein et al. 2008). The clusters 
studied in the present work have distances in the range of ~ 0.6 to ~ 3 kpc. The polarimetric 
results obtained can be used to study the properties of the dust towards the anticenter 
direction {I ~ 160° - 175°) of the Galaxy. Based upon the earlier discussion using Q v — Uy 
and (U—B)/ (B — V) diagrams, in Fig. [5] we show a combined Uy — Qy diagram for all the five 
clusters, namely, NGC 2281 (inverted triangles), NGC 1664 (squares), NGC 1960 (triangles), 
Stock 8 (circles) and NGC 1893 (star symbols). Open and filled symbols are used to represent 
members and non-members based upon our analysis. The dust free environment of the solar 
neighbourhood (shown with rectangular box) is represented by Qy = and Uy = 0, while 
any other point on this plot represents the direction of the polarization vector seen towards 
that direction from the Sun. Fig. [5] suggests that the degree of polarization of stars is found 
to increase with the distance to the clusters, as also noticed in Fig. |2j This is consistent with 
the fact that the degree of polarization increases with the column density of dust grains 
lying in front of the stars that are relatively well aligned. The close proximity of the points 
of NGC 1664 and NGC 1960 in Uy — Qy plane is consistent with the fact that they are 
located approximately at similar distances (1.2 and 1.3 kpc, respectively). 
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Figure 5. The Uv versus Qv diagram for stars lying in the cluster regions studied, namely, NGC 2281 (inverted triangles), 
NGC 1664 (squares), NGC 1960 (triangles), Stock 8 (circles) and NGC 1893 (star symbols). The filled symbols represent field 
stars in the respective cluster regions. The rectangular box covered by Qy = and Vy = is the dust free environment of the 
solar neighbourhood. 



4.2 The Serkowski law 

The wavelength dependence of polarization towards many Galactic directions follows the 
Serkowski's law of interstellar polarization (Serkowski 1973; Coyne, Gehrels & Serkowski 
1974; Wilking, Lebofsky, & Rieke 1982): 

P\ = P m a X exp[-K \n 2 (X max /X)} (3) 

Where P\ is the percentage polarization at wavelength A and P max is the peak polarization, 
occurring at wavelength \ ma x- The \ ma x is a function of optical properties and characteristic 
particle size distribution of aligned grains (Serkowski, Mathewson & Ford, 1975; McMillan, 
1978). The value of P max is determined by the column density, the chemical composition, size, 
shape, and alignment efficiency of the dust grains. The parameter K, an inverse measure of 
the width of the polarization curve, was treated as a constant by Serkowski et al. (1975), who 
adopted a value of 1.15 for all the stars. The Serkowski's relation with K—1.15 provides an 
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adequate representation of the observations of interstellar polarization between wavelengths 
0.36 and 1.0 /im. If the polarization is produced by aligned interstellar dust grains, the 
observed data will follow equation (ED and hence we can estimate P ma x and X ma x for each 
star. The P max and X max are obtained using the weighted least-squares fitting to the measured 
polarization in BV(RI) C bands to equation (j3J) by adopting K—1.15. We have also computed 
the parameter crj^l (the unit weight error of the fit) for each star that quantifies the departure 
of the data points from the standard Serkowski's law. Because of the weighting scheme, the 
values of <j\ should not exceed 1.5, but if they do, it implies that the stars have intrinsic 
polarization (e.g., Feinstein et al. 2008, Medhi et al 2007, 2008, 2010). The \ max values can 
also be used to infer the origin of the polarization. The stars with \ max much lower than the 
average value of the interstellar medium (0.545 /xm; Serkowski, Mathewson & Ford 1975) 
may have an intrinsic component of polarization. Another parameter to infer the presence 
of intrinsic polarization or polarization angle rotation along the line of sight (Coyne 1974, 
Martin 1974), is the dispersion of the polarization angle for each star normalized by the 
average of the polarization angle errors, e (Marraco et al. 1993; Orsatti et al. 2007; Feinstein 
et al. 2008). 

The estimated values of P ma x, X max , &i and e for the stars towards NGC 1893 are given 
in Table [9j The star identifications are same as given in Table HJ The weighted mean values 
of the P max and A max are found to be 2.59 ± 0.02% and 0.55 ± 0.01 /im respectively. The 
estimated X ma x is quite similar to the value corresponding to the general interstellar medium 
(0.545 /an, Serkowski et al. 1975). Using the relation Ry = 5.6 x X max (Whittet & Van Breda 
1978), the value of Ry, the total-to-selective extinction, comes out to be 3.08±0.05, which is 
in agreement with the average value {Ry = 3.1) for the Milky Way Galaxy, indicating that 
the size of the dust grains within the cluster NGC 1893 is normal. Similar conclusion was 
drawn by Sharma et al. (2007) using (A — V)/(B— V) two colour diagrams. In Fig. [6] we show 
(Ji versus P max (upper panel) and e versus \ max (lower panel) plots. The criteria mentioned 
above indicate that majority of the stars do not show evidence of intrinsic polarization. 
However a few stars #72, 149, 199, 209 and 215 show deviation from the general distribution. 

Star #199 is an Ha emission star and its location in the V/(V — I) colour-magnitude 
diagram (CMD) (Fig. [7]) reveals that it should be a pre-main sequence (PMS) star. The 

The values of a\ for each star are computed using the expression cr^ — ^ ^ (V \ j 'tpx^) 2 j '{rn — 2); where yyi is the nnmber of 
colours and r x = P\ — Pmax exp[— K In 2 ''(Xmax / 'A) . 
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Figure 6. Top panel is the plot between a\ and P m ax, whereas the bottom panel is the plot between e versus \ m ax- 

location of star #209 in the (U — B)/(B — V) colour-colour diagram (Fig. HJ) suggests that 
it could be a foreground field star. The V/ (V — I) CMD reveals that star #215 should be a 
PMS star. Star #123 and 199 are identified as Ha emission stars S1R2N44 and S1R2N56, 
respectively (Neguerela et al. 2007). Stars #199 and 215 are likely to be PMS stars and ac- 
creting matter from their circumstellar material. The accreted material, probably distributed 
in an asymmetric disk geometry, might be responsible for the observed intrinsic component 
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Figure 7. Upper panel is the V versus (V — I) colour-magnitude diagram for the NGC 1893 region. The continuous, dotted 
and dashed curves represent the ZAMS from Girardi et al. (2002) and PMS isochrones for 1 and 5 Myr by Siess et al. (2000) 
respectively. The isochrones are adjusted for a distance of 3.2 kpc and E(B — V)=0.40 mag. The stars having Ha emission are 
shown by crosses. The lower panel is the enlarged view of the colour-magnitude diagram shown in the inset of the upper panel. 

of polarization. The computed values of e show no significant rotation (except for star #149) 
in polarization angles reiterating the fact that the dust component responsible for the mea- 
sured polarization is well aligned. Star #149 is a highly reddened MS star (see. Figs EH [7] 
and [8]) as it is embedded at the eastern edge of the wind blown bubble (cf. Sharma et al. 
2007). 

4.3 Polarization efficiency 

The degree of polarization produced for a given amount of extinction (or reddening) is re- 
ferred as the polarization efficiency of the intervening dust grains. The polarization efficiency 
depends mainly on the orientation of the magnetic field along the line of sight, the mag- 
netic field strength and the degree of alignment of the dust grains. Mie calculations place 
a theoretical upper limit of P ma x = 43 x E(B — V) (Whittet 1992 and references there 
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Figure 8. (J — H) versus (H — K) colour-colour diagram for the NGC 1893 region. The data is taken from the Two-Micron 
All-Sky Survey (2MASS) Point Source Catalog (PSC) (Cutri ct al. 2003). 2MASS data has been converted in to California 
Institute of Technology (CIT) system using the relations provided by Carpenter (2001). The theoretical tracks for dwarfs and 
giants arc drawn (Bessell & Brett, 1988). Reddening vectors are also drawn (Cohen et al. 1981). The location of Be stars (cf. 
Dougherty et al. 1994), and the location of Herbig Ae/Be stars (cf. Hernandez ct al. 2005) are also shown. The stars having 
Ha emission are shown by crosses. 

in) on the polarization efficiency by an infinite cylinder with the diameter comparable to 
the wavelength of the incident light and with their long axes parallel to each other and 
perpendicular to the line of sight. The empirical upper limit relation for the polarization 
efficiency resulting from the studies of reddened Galactic stars, assuming normal interstellar 
material characterized by Ry = 3.1 (Serkowski, Mathewson & Ford 1975) is found to be 
Pmax = 9 x E(B — V) . The fact that the maximum observed polarization efficiency is found 
to be less by a factor of 4.8 than that expected from theory, implies that the alignment of 
dust grains is not perfect. The reason could be the presence of various components of mag- 
netic fields oriented differently along the line of sight and/or the grains are only moderately 
elongated rather than infinite cylinder. 

The UBV photometric data by Sharma et al. (2007) and Massey et al. (1995) have been 



20 C. Eswaraiah et al. 




0.0 0.2 0.4 0.6 0.8 1.0 

E(B-V) 

Figure 9. The polarization efficiency diagram. Filled and open circles are the stars distributed in the northern and southern 
regions of the cluster respectively. The solid line represents the empirical upper limit relation for the polarization efficiency of 
Pmax = 9 X E(B — V) (Serkowski et al. 1975). The dashed line represents the relation P m ax = 5 X E(B — V) (Serkowski et al. 
1975) and the dashed-dotted line represents the relation P max = 3.5 X E(B - V) - 8 by Fosalba et al. (2002). 

used to estimate the reddening E(B — V). The reddening of individual stars having spectral 
type earlier than AO, has been derived using the Q method (Johnson & Morgan 1953) and 
the values are given in Table |9j The (J — H)/(H — K) near infrared (NIR) colour-colour 
diagram (Fig. [8]) reveals that stars #48, 72, 197, 310, 355 and 356 have spectral type later 
than AO. The (U — B)/(B — V) colour-colour diagram (Fig. H]) indicates that these stars 
follow ZAMS reddened by E(B — V)=0.2 mag, hence we assign reddening of E(B — V)=0.20 
mag. 

In Fig. [91 we present P ma x versus E(B — V) for stars towards NGC 1893 region. The 
continuous line shows the empirical upper limit for the polarization efficiency given by 
Pmax = 9 x EiB — V) (Serkowski et al. 1975), assuming normal interstellar material charac- 
terized by Ry = 3.1. The recent estimate of the average efficiency by Fosalba et al. (2002), 
which is valid for E(B — V)< 1.0 mag, is represented by the relation P m ax = 3.5 x E{B — V) 8 
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Figure 10. The radial variation of P m ax (upper panel) and E(B — V) (lower panel) for the cluster members and the stars (# 
108, 123, 196, 199 and 213) with uncertainty in their membership. Filled circles are the stars distributed in the northern part, 
whereas open circles are in the southern part of the cluster NGC 1893. Most of the stars falling within the circle at a distance 
of 5.'5 are distributed nearer the two nebulae Sim 129 and 130. 



and is shown by the dash-dotted line. For comparison, the average polarization efficiency 
relation, P ma x = 5 x E(B — V) (Serkowski et al. 1975), is also drawn using the dashed line. 
A majority of the points are found to lie between the continuous line and above the dashed 
line, indicating that the dust grains towards NGC 1893 have higher polarization efficiency 
in comparison to the the general ISM. 
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Figure 11. Cumulative distribution of Pmax for the stars, with E(B — V) values between 0.4-0.6 mag, distributed in northern 
and southern parts of the cluster NGC 1893. The lines connecting filled and open circles are the cumulative distributions for 
northern and southern stars respectively. 

4.4 Spatial variation of P max and E(B — V) 

The upper panel of Fig. [TOl shows P max for the cluster members and the stars with uncertainty 
in their membership 108, 123, 196, 199 and 213), as a function of the radial distance from 
the ionization source (see Fig.[TJ HD 242935; shown by the plus symbol). Stars located in the 
northern and southern parts of the cluster are shown by filled and open circles respectively. 
The distribution of P ma x within a radial distance of ~ 5 from the ionization source shows a 
decreasing trend with the increase in radial distance. The distribution of P ma x shown in Fig. 
[TOl also reveal that P ma x m the southern region is relatively higher than that of the northern 
region. There is a significant increase in the P max values at ~ 5f5. The increase in P max is 
mainly due to the YSOs distributed around the nebulae Sim 129 and 130, located towards 
the North-East direction of HD 242935 at a radial distance of ~ 6 . 

The lower panel of Fig. [10] shows radial variation of E(B — V). The distribution of 
E(B — V) reveals a clumpy nature of the gas/dust distribution. The polarization efficiency 
in the northern region of the cluster at a radial distance ~ O'.b, assuming a normal value of 
i?y=3.1, comes out to be ~ 6.0, and decreases to ~ 5.0 at the boundary (~ 4 ) of the cluster. 
The ratio P max /E(B — V) depends mainly on the alignment efficiency, magnetic strength 
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and amount of depolarization due to radiation traversing through more than one medium. 
The polarization efficiency P\/A\ also depends on the particle shape (Voshchinnikov & Das 
2008). 

The weighted mean of P ma x for stars having 0.40 ^ E(B — V) ^ 0.60 lying in the northern 
and southern regions of the cluster is estimated to be as 2.64±0.04% and 2.77±0.07%. 
Fig. HH shows the cumulative distribution of P ma x in both regions, for stars having 0.40 
^ E(B — V) ^ 0.60 mag, which indicates that P ma x values in the southern region are 
systematically higher than those in the northern region. Kolmogorov-Smirnov test indicates 
that these distributions are different at 80% confidence level. 

The weighted mean values of polarization angle (8y) for northern and southern regions 
are 161 ±4° and 162 ±6°, which indicate that there is no difference in 9y in the two regions. 

5 DUST COMPONENTS RESPONSIBLE FOR THE OBSERVED 
POLARIZATION 

The degree of polarization of a star increases as a function of distance due to the presence 
of a column of aligned dust grains along the pencil beam of radiation from the star. The 
degree of polarization shows a sudden jump if the radiation from the stars encounter a dust 
layer located at a certain distance between the star and the observer. The number of such 
sudden jumps is characterized by the number of dust layers encountered by the radiation 
along its path and the relative magnetic field orientations in the dust layers. 

To understand the distribution of dust layers towards the direction of NGC 1893, we 
selected stars from a region of 10° radius around the cluster with both polarization measure- 
ments (Heiles 2000) and Hipparcos parallaxes (van Leeuwen, 2007) available. In Fig. [121 we 
show the degree of polarization versus distance (upper panel) and polarization angle versus 
distance (lower panel) plots for the selected stars (filled circles). Stars observed towards the 
clusters NGC 2281, NGC 1664, NGC 1960, Stock 8 and NGC 1893 are also shown using 
inverted triangles, squares, triangles, open circles and star symbols respectively. Two sig- 
nificant jumps in the values of Py, one at ~ 170 pc and another at ~ 360 pc, are clearly 
evident in the Py versus distance plot (upper panel). 

For distances < 170 pc, the 9y values show a scattered distribution between ~ 50° — 130°. 
At ~ 170 pc, the 9y values show a dip to ~ 20°. A sharp rise in 8y values from ~20° to 160° 
is evident at a distance of ~ 360 pc. The stars located beyond ~ 1000 pc show an average 
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polarization angle of 9y ~ 163°. On the basis of Fig. [121 we infer that a dust layer at ~ 170 
pc contributes ~ 0.3 — 0.9% to the polarization and the dust grains are aligned towards 9y 
~ 20°. Another dust layer with 9y ~ 160 — 170° at a distance ~ 360 pc further contributes 
~ 0.3 — 1.3% (i.e., total ~ 1.2 to ~ 2.2 %). The clusters located at distance > 1000 pc show 
Py in the range of ~ 1 — 4%. The Py > 2.2% in these clusters must be due to intracluster 
medium. The optical properties derived in Sections 14.21 and 14.31 should be a combined effect 
of dust layers at d ~ 170 pc, ~ 360 pc and intracluster medium. 

Neckel and Klare (1980) have studied Ay distribution in the Galactic plane with |6| ^ 7.6° 
using extinction and distances computed for individual stars. The Ay map towards the 
direction of NGC 1893 by Neckel and Klare (1980) reveals that Ay increases with distance 
up to ~ 2 kpc. Beyond ~ 2 kpc Ay is found to be rather constant, indicating that the dust 
contribution in the distance range 2 — 3 kpc is negligible. Here, it is worthwhile to mention 
that the mean polarization Py in the Stock 8 (2.4±0.6%) and NGC 1893 (2.6±0.7%) is 
almost the same, which further confirms that the contribution of dust in the distance range 
~ 2-3 kpc is negligible. Further, the observed mean polarization values are consistent with 
the E(B — V) values mentioned in Tabled] 

6 CONCLUSIONS 

We have made polarimetric observations of the open cluster NGC 1893 in B, V, R c and 
I c bands. Three additional open clusters towards the direction of NGC 1893, namely, NGC 
2281, NGC 1960 and Stock 8 along with NGC 1664 toward anticenter direction of the Galaxy 
were observed in V^-band only. The main aim of the study is to investigate the properties 
of dust grains towards the anticenter direction of the Galaxy between I ~ 160° to ~ 175°, 
using stars of open clusters located in the distance range of 600 pc (NGC 2281) to 3.2 kpc 
(NGC 1893). 

The stars located at distances ^ 170 pc show a large scatter in 6y values with a polar- 
ization of ~ 0.1%. The degree of polarization is found to increase with the distance. The 
distribution of Py and 9y as a function of distance reveals two dust layers at ~ 170 pc and 
~ 360 pc. The first dust layer is characterized by its polarized components as Py ~ 0.3—0.9% 
and 9y ~ 20 — 50°. Both dust layers produce a maximum combined polarization up to Py 
~ 2.2%. Polarization values higher than the ~ 2.2% should be due to intracluster medium. 

The magnetic field orientation remains unchanged within the Galactic plane |6| < 2°. 
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Figure 12. The change in Py(%) and 9v{°) as a function of distance (pc) is plotted in the upper and lower panels for stars 
(filled circles) selected from a circular region of 10° around NGC 1893. The P% and 8° values are obtained from the Heiles 
(2000) catalog and the distances are obtained from Hipparcos parallaxes (van Leeuwen, 2007). The vertical solid lines are drawn 
at 170 pc and 360 pc. The results from our study of the five clusters are also shown. The clusters are represented by the same 
symbols as in Fig. [5] 



The estimated mean value of the polarization angle (assuming a Gaussian distribution) of 
the four clusters (NGC 1664, NGC 1960, Stock 8 and NGC 1893) comes out to be ~ 163° 
with a standard deviation of 6°. This small dispersion in polarization angle could be due to 
the presence of a uniform dust layer beyond 1 kpc. Present observations reveal that in case 
of NGC 1893, the foreground two dust layers, in addition to the intracluster medium, seems 
to be responsible for the polarization effects. 

The estimated mean values of Py for the two clusters namely, Stock 8 and NGC 1893, 
imply that the ISM located between 2-3 kpc distance has negligible contribution towards 
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extinction as well as the observed polarization. Present polarimetric results are consistent 
with the reddening distribution given by Neckel and Klare (1980). 

The weighted mean of the P ma x and \ max values for NGC 1893 are found to be 2.59 ± 
0.02% and 0.55 ± 0.01 /iin respectively. The estimated \ max is quite similar to that of the 
general interstellar medium. The value of R v using the relation R v = 5.6 x \ max (Whittet & 
Van Breda 1978) is found to be close to the average value of 3.1 for the Milky Way Galaxy, 
implying that the average size of the dust grains within the cluster NGC 1893 is similar 
to the general interstellar medium. We also identified four candidate stars with intrinsic 
polarization in NGC 1893. 

The radial distribution of P ma x within the cluster shows a decrease of P max towards 
the outer region of the cluster. The P ma x and the polarization efficiency is found to be 
higher towards the southern region of the cluster. We have shown that the polarization 
measurements in combination with (U — B)/(B — V) colour-colour diagram provide a good 
tool to determine the membership in a cluster. 
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Table 3. Observed polarized and unpolarizcd standard stars 





P±e P (%) 


e±e e (°) 


P±e P (%) 


e±e e n 




Our work 


Schmidt et al. (2002) 


Polarized standard stars 


2008 November 8 & 9 
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B 


6.32± 0.22 


97.9±1.0 


6.34 ±0.04 


98.14 ±0.16 


V 


6.86± 0.16 


97.3±0.7 


6.70 ±0.01 


98.09 ±0.07 


Rc 


6.43± 0.13 


96.0±0.6 


6.43 ±0.02 


98.14±0.10 


la 


5.86± 0.11 


97.0±0.5 


5.80 ±0.02 


98.26±0.11 






BD+64° 106 




B 


5.73± 0.21 


97.3±1.0 


5.51 ±0.09 


97.15±0.47 


V 


5.68± 0.19 


96.7±1.0 


5.69 ±0.04 


96.63±0.18 


Rc 


5.12± 0.17 


97.4±0.9 


5.15 ±0.10 


96.74±0.54 


Ic 


4.66± 0.21 


97.8±1.3 


4.70 ±0.05 


96.89±0.32 
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B 


5.65± 0.22 


97.4±1.1 


5.51 ±0.09 


97.15±0.47 
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5.49± 0.20 


96.9±1.0 


5.69 ±0.04 


96.63±0.18 


Rc 


5.33± 0.18 
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4.70 ±0.05 


96.89±0.32 






HD204827 




B 
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Table 4. Observed B, V , (RI) C polarization and polarization angles for stars towards NGC 1893 
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2.8 


2.65 


± 


0.36 


165.8± 


3.7 


2 


.38 


± 


0.46 


160.6 


± 


5.3 


327 


11.62" 


1 


.25 


± 


0.22 


161.7 


± 


4.7 


1.01 


± 


0.21 


164.3 


± 


5.3 


1.06 


± 


0.26 


156.0± 


6.5 


1 


.05 


± 


0.34 


168.4 


± 


8.5 


329 


11.95" 


1 


.20 


± 


0.26 


152.7 


± 


5.7 


0.84 


± 


0.24 


160.5 


± 


7.4 


1.03 


± 


0.32 


148.4± 


8.0 


1 


.49 


± 


0.42 


164.5 


± 


7.6 


355 


13.46" 


1 


.61 


± 


0.53 


154.8 


± 


8.8 


1.94 


± 


0.47 


168.2 


± 


6.6 


2.01 


± 


0.53 


168.7± 


7.2 


2. 


.01 


± 


0.64 


170.8 


± 


8.6 


356 


14.35" 


2. 


.57 


± 


0.80 


160.7 


± 


8.6 


2.01 


± 


0.67 


164.0 


± 


8.7 


3.23 


± 


0.58 


157.3± 


5.0 


1 


.85 


± 


0.60 


150.6 


± 


8.7 



t: Cuffey and Shapley (1937) 
*: Sharma ct al. (2007) 
a : Massey et al. (1995) 
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Table 5. Polarimctric results for 14 stars observed towards NGC 2281. 



1A t 

Id. 1 


IO \ 

a{ ) 


<H ) 


V 


(B - 


- V)* 


It t 

(U - 


- B )* 




Pv 




0v 




i\/rttt 
M p ' 




(J 2000) 


(J2000) 


(mag) 












iw \ 

(%) 




( ) 






(!) 


( 2 ) 


(3) 


( 4 ) 


( 


5 ) 


( 


6) 




(J) 




(8) 


(9) 


(10) 


49 


102.03038 


41.077494 


10.56 


0. 


24 





.16 


1 


.07± 


.17 


16.16 ± 4.28 


M 


M (99) 


53 


102.06064 


41.112788 


11.02 


0. 


33 


0. 


.13 





,91± 


22 


19.52 ± 6.14 


M 


M (99) 


55 


102.06337 


41.073019 


8.89 


0. 


98 





69 


1 


,00± 


08 


13.77 ± 2.10 


NM 


M (99) 


56 


102.06596 


41.095871 


11.43 


1. 


24 


1. 


19 


1 


,04± 


25 


6.23 ± 6.39 


NM 


NM (00) 


57 


102.07032 


41.130267 


10.61 


0. 


23 





15 





90± 


18 


6.39 ± 5.05 


M 


M (99) 


58 


102.07208 


41.083118 


9.45 


0. 


12 





06 





.73± 0. 


10 


9.65 ± 3.53 


M 


M (99) 


60 


102.08048 


41.041193 


12.71 


0. 


56 


0. 


07 


1 


.41± 0. 


.45 


12.29 ± 8.47 


M 


M (99) 


71 


102.11072 


41.024925 


11.16 


0. 


53 


0. 


05 





.78± 0. 


23 


24.01 ± 7.43 


M 


NM (00) 


74 


102.1166 


41.081101 


9.54 


1. 


10 


1. 


03 





.56± 0. 


11 


22.08 ± 4.62 


NM 


NM (00) 


78 


102.13344 


41.108674 


10.63 


0. 


22 





.15 





60± 


.18 


23.54 ± 7.22 


M 


M (99) 


79 


102.14759 


41.090775 


10.30 


0. 


17 





.11 





.67± 


15 


12.83 ± 5.60 


M 


M (99) 


82 


102.15634 


41.10563 


10.08 


0. 


14 





10 





.85± 


14 


18.06 ± 4.11 


M 


M (99) 


86 


102.17038 


41.066299 


8.62 


0. 


22 


0. 


15 





.93± 0. 


07 


19.43 ± 1.92 


M 


M (99) 


91 


102.18262 


41.086379 


11.41 


0. 


35 





12 





.94± 0. 


25 


25.04 ± 6.97 


M 


M (99) 



t : Vasilevskis & Balz (1959) 

*: Pesch (1961) except star no # 60 for which V magnitude, (B — V) and (U — B) colours are taken from Yoshizawa (1978) 

tt; Present Work 

ttt : Vasilevskis & Balz (1959) 



Polarimetric study towards NGC 1893 



Table 6. Polarimetric results on 


27 stars observed towards NGC 1664. 












Id.t 


a(°) 


S(°) 


V 


(B-V) 


(U-B) 


Pv 


e v 




M„ tt 


M p ttt 




(J2000) 


(J2000) 


(mag) 






(%) 


(°) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


(10) 


1 


72 783342 


43 6651 61 


13.50™ 


1.61° 


1.68° 


2 44 + 35 

Zi . *-±^t _1_ U.OJ 


165 


± 


4 


NM 




4 


79 7fiQ^nn 


43 671 694 


11.75^ 


0.32° 


0.25° 


1.41 _t 0.15 


178 


± 


3 


M 


JVL {00 ) 


7 


72 824534 


43 671 Q85 


12.25° 


0.30° 


0.24° 


1.14 ± 0.18 


165 


± 


4 


M 


T\/T fQC\\ 

JVL [ov ) 


9 


72 81 4747 


43 6801 35 








O 7c: _|_ n 1 a 


164 


± 


1 


NM 




12 


72 749866 


43 684928 


13.01" 


0.24° 


0.31° 


1 0^ + 27 

.L.ZitJ _1_ \J.Zj\ 


170 


± 


6 


M 


JVL [\)2 ) 


13 


79 7608Q5 
i z, . i uuoyu 


43 676394 


14.36° 


1.09° 


0.71° 


O KQ _|_ f| 4Q 

Z1.0& _tz u.iy 


163 


± 


5 


IN M 




14 


72 750261 


43 676335 


12.79™ 


0.51° 


0.28° 


1 25 + 26 


170 


± 


5 


M 


JVL (ylj 


21 


72 771228 

i Z, ■ i 1 liiU 


A% 648455 


11.816 


0.29° 


0.26° 


1 ,39 + 15 


173 


± 


Q 
O 


"v t 

M 


M [p2 ) 


22 


72 779262 

I Zj . ) I C7ZjVJZj 


43 640397 


14.00™ 


0.54° 


0.38° 


2 S3 + 42 

Zj . cjo _i_ \j.'-tZj 


171 


± 


5 


IN M 


TV/T 

M {00) 


27 


79 81 6757 


43 6569Q1 


13.87° 


0.51° 


0.21° 


9 1 -I- 3Q 
^ . ±u _tz u.oy 


175 


± 


5 


JNIM 


INI M {00 ) 


35 


72 804006 


43 700468 


12.50 c 


0.37° 


0.25° 


1 15 + 021 


170 


± 


5 


M 


M (ylj 


37 


79 7Q0Q09 


43 708346 


10.93° 


0.34° 


0.34° 


1 9zL -1- 1 

_L . Zj^i _TZ U.1U 


174 


± 


2 


M 


M (90) 


47 


72 754574 


43 703638 








3 18 H - 54 


160 


± 


5 


IN M 




55 


79 7941 55 

i Zi. i Zi L ±\-0<) 


43 669579 


11.06° 


1.23° 


1.00° 


1.71 i 0.11 


169 


± 


2 


M 


M (71) 


56 


72.737611 


43.657761 


12.66° 


0.35° 


0.26° 


1.40 ± 0.23 


177 


± 


A 

4 


M 


TV. \ /m A 

M (yij 


67 


72.857380 


43.672892 








1.21 ± 0.38 


171 


± 


8 


Ad 


M (9z!J 


68 


72.849878 


43.675125 








1.57 ± 0.27 


168 


± 


5 


M 


M (87) 


75 


72.832061 


43.704486 


11.25° 


1.06° 


0.83° 


1.47 ± 0.13 


173 


± 


2 


M 


M (92) 


76 


72.823740 


43.720916 








1.53 ± 0.36 


177 


± 


6 


M 


M (87) 


77 


72.814271 


43.709793 


12.84 c 


0.49° 


0.25° 


1.40 ± 0.25 


174 


± 


5 


M 


M (65) 


79 


72.765348 


43.727229 


13.75° 


0.54° 


0.50° 


2.06 ± 0.38 


165 


± 


5 


NM 




87 


72.700196 


43.654016 








4.09 ± 0.80 


171 


± 


5 


NM 




90 


72.712849 


43.643313 


13.03° 


0.51° 


0.30° 


1.47 ± 0.29 


170 


± 


5 


M 


NM (17) 


113 


72.735822 


43.721862 


12.97° 


0.53° 


0.27° 


1.45 ± 0.29 


175 


± 


5 


M 


M (84) 


154 


72.780622 


43.746898 








1.26 ± 0.31 


172 


± 


6 


M 




158 


72.719637 


43.739690 


14.79° 


0.73° 


0.13° 


1.78 ± 0.50 


166 


± 


7 


M 




164 


72.691065 


43.720635 








0.97 ± 0.27 


26 ± 7 


NM 


NM (00) 



t; Larsson-Leander (1957) 

°: Hoag ct al. (1961) Publ. Us. Nav. Obs. XVII part VII, 347 
b : Purgathofer (1964) 
c : Hoag et al. (1961) 
tt: Present Work 
ttt ; Dias ct al. (2006) 
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Table 7. Polarimctric results on 15 stars observed towards NGC 1960. 



TJ t 

Id. 1 


a{ ) 


o( ) 


V 


(B - V)* 


(U -B)* 




Pv 


0v 




\S tt 


i/i ttt 
M p 11 1 




(J 2000) 


(J2000) 


(mag) 








(%) 


(°) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 




(9) 


(10) 


13 


84.063492 


34.120250 


10.78 


0.12 


-0.26 


1 


.12 ± 0.17 


160 


± 


4 


M 


M (94) 


16 


84.065508 


34.143746 


8.86° 


-0.00 a 


-0.66 a 


1 


.17 ± 0.07 


159 


± 


2 


M 


M (94) 


17 


84.060257 


34.139917 


12.40 


0.25 


0.18 


1 


.16 ± 0.35 


162 


± 


8 


M 


M (94) 


23 


84.095740 


34.175955 


8.96 a 


0.01 a 


-0.68 a 


1 


.12 ± 0.07 


162 


± 


2 


M 


M (94) 


33 


84.118160 


34.122944 


11.85 


0.12 


-0.07 


1 


.15 ± 0.27 


155 


± 


6 


M 


M (94) 


38 


84.099176 


34.099402 


9.92" 


0.05 a 


-0.49 a 


1 


.17 ± 0.11 


156 


± 


2 


M 


M (94) 


41 


84.083418 


34.103419 


12.37 


0.19 


0.14 


1 


.28 ± 0.35 


160 


± 


7 


M 


M (94) 


44 


84.047893 


34.118557 


11.35 


0.09 


-0.36 


1 


.38 ± 0.22 


163 


± 


4 


M 


M (94) 


55 


84.081130 


34.204218 


11.61 


0.10 


-0.31 


1 


.44 ± 0.24 


147 


± 


5 


NM 


M (94) 


56 


84.087224 


34.211974 


12.39 


0.46 


0.35 


1 


.15 ± 0.26 


157 


± 


6 


M 


M (93) 


61 


84.132976 


34.179480 


9.14" 


0.01° 


-0.66 a 





.81 ± 0.08 


155 


± 


2 


NM 


M (94) 


77 


84.069337 


34.083562 


12.10 


0.17 


0.06 


1 


.50 ± 0.31 


162 


± 


6 


NM 


M (94) 


87 


83.996719 


34.174535 


10.62 


0.07 


0.95 


1 


.34 ± 0.16 


161 


± 


3 


M 


M (94) 


91 


84.040637 


34.193293 


10.34 


0.01 


-0.50 


1 


.25 ± 0.14 


165 


± 


3 


M 


M (94) 


92 


84.024453 


34.201113 


10.93 


0.03 


-0.49 


1 


.16 ± 0.18 


164 


± 


4 


M 


M (94) 



t : Boden (1951) 

*: Sharma et al. (2006) 

a : Johnson H.L., Morgan W.W. (1953) 

tt; Present Work 

ttt ; Dias et al. (2006) 



Table 8. Polarimctric results on 21 stars observed towards Stock 8 



Id.t 


«(°) 


6(°) 


V 


(B-V)(*) 


(U-B)(*) 


Pv 


0v 




Mptt 


Mpttt 




(J2000) 


(J2000) 


(mag)* 






(%) 


(°) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


(10) 


11 


82.004547 


34.450581 


11.50 


0.13 


0.17 


1.57 ±0.14 


163 ± 


2 


NM 


M (66) 


12 


82.009061 


34.404500 


11.76 


-0.36 


0.28 


2.41 ±0.15 


165 ± 


2 


M 


M (66) 


13 


82.018190 


34.489546 


11.71 


-0.47 


0.30 


2.45 ±0.15 


167 ± 


2 


M 


M (64) 


16 


82.038176 


34.473935 


11.89 


-0.40 


0.29 


2.41 ±0.16 


151 ± 


2 


M 


M (65) 


19 


82.058027 


34.438897 


13.01 


1.22 


1.40 


2.17 ±0.28 


157 ± 


4 


NM 


M (77) 


141 


82.082569 


34.420441 


12.19 


-0.24 


0.51 


2.60 ±0.19 


165 ± 


2 


M 


M (57) 


160 


82.002940 


34.493490 


12.63 


0.15 


0.60 


1.62 ±0.23 


163 ± 


4 


NM 


NM (02) 


173 


82.082563 


34.455984 


12.81 


0.37 


0.31 


1.60 ±0.25 


159 ± 


4 


NM 


NM (47) 


181 


82.064303 


34.423842 


12.87 


-0.32 


0.33 


2.39 ±0.25 


158 ± 


3 


M 


M (68) 


208 


82.040118 


34.446178 


13.18 


-0.14 


0.30 


2.62 ±0.28 


162 ± 


3 


M 


M (74) 


211 


81.996346 


34.444053 


13.21 


0.17 


0.28 


1.93 ±0.30 


159 ± 


4 


NM 


NM (19) 


218 


82.004829 


34.385705 


13.29 


-0.02 


0.34 


2.61 ±0.31 


162 ± 


3 


M 


NM (43) 


231 


82.028692 


34.394823 


13.36 


-0.14 


0.31 


3.03 ±0.33 


157 ± 


3 


M 


M (72) 


258 


81.969263 


34.389477 


13.53 


0.24 


0.68 


2.57 ±0.35 


161 ± 


4 


M 


NM (00) 


275 


82.090470 


34.409075 


13.69 


-0.17 


0.62 


2.37 ±0.39 


136 ± 


5 


? 




294 


82.003915 


34.403745 


13.81 


0.38 


0.46 


2.72 ±0.40 


163 ± 


4 


M 


M (52) 


320 


81.993454 


34.470836 


13.93 


0.27 


0.62 


2.20 ±0.45 


156 ± 


6 


M 


NM (01) 


325 


81.962444 


34.428255 


13.96 


-0.08 


0.82 


4.18 ±0.44 


165 ± 


3 


? 


M (71) 


333 


82.066055 


34.399693 


14.00 


0.32 


0.59 


2.40 ±0.45 


163 ± 


5 


M 




365 


82.001149 


34.457243 


14.18 


0.37 


0.41 


2.17 ±0.46 


165 ± 


6 


NM 




425 


81.991192 


34.498474 


14.39 


0.19 


0.55 


2.36 ±0.61 


158 ± 


7 


M 


NM (04) 



t: Mayer (1964) 
*: Jose et al. (2008) 
tt: Present Work 
ttt; Dias et al. (2006) 



Polarimetric study towards NGC 

Table 9. The P ma x, Xmax, <ri & e for the observed 44 stars towards NGC 1893 



Id.t 


E(B - V)* 


a (°) 


5(°) 


p 


ax 


± 


e 




ax 


± 


e 




e 


Mptt 


Mpttt 




(mag) 


(J2000) 


(J2000) 




(%) 




(jim) 












(1) 


(2) 


(3) 


(4) 






(5) 








(6) 






(7) 


(8) 


(9) 


(10) 


1 


0.53 


80.683340 


33.440739 


2 


50 


± 





08 


0.53 


± 





04 


0.74 


1.31 


M 


M (78) 


2 


0.54 


80.686558 


33.443501 


3 


11 


± 





13 


0.54 


± 





05 


0.69 


1.23 


M 


M (79) 


3 


0.45 


80.706563 


33.447926 


2 


41 


± 





08 


0.58 


± 





04 


0.91 


1.38 


M 


M (78) 


5 


0.46 


80.725969 


33.445042 


2 


52 


± 





26 


0.52 


± 





10 


0.59 


0.47 


M 


NM (9) 


9 


0.52 


80.706864 


33.425991 


2 


76 


± 





18 


0.53 


± 





07 


0.89 


1.25 


M 


M (75) 


10 


0.46 


80.705199 


33.428371 


2 


61 


± 





28 


0.70 


± 





14 


1.03 


0.41 


M 


NM (0) 


13 


0.52 


80.692334 


33.422447 


2 


80 


± 





22 


0.62 


± 





10 


0.54 


0.54 


M 


M (79) 


15 


0.48 


80.696872 


33.418720 


2 


89 


± 





25 


0.57 


± 





10 


0.83 


0.92 


M 


M (61) 


33 


0.50 


80.688066 


33.406540 


2 


32 


± 





20 


0.56 


± 





10 


0.62 


1.27 


M 


M (71) 


48 


0.20 a 


80.699353 


33.476135 


1 


76 


± 





14 


0.51 


± 





06 


1.46 


0.17 


NM 


M (60) 


59 


0.47 


80.741815 


33.443459 


2 


26 


± 





13 


0.60 


± 





07 


0.72 


0.98 


M 


M (79) 


62 


0.47 


80.739711 


33.433380 


2 


41 


± 





16 


0.56 


± 





07 


0.82 


0.65 


M 


M (74) 


72 


0.20 a 


80.699203 


33.368629 


1 


86 


± 





40 


0.87 


± 





23 


0.82 


1.05 


NM 


NM (27) 


73 


0.27 


80.688608 


33.371365 


1 


80 


± 





11 


0.53 


± 





06 


0.58 


0.50 


NM 


M (79) 


82 


0.70 


80.654111 


33.386971 


3 


32 


± 





26 


0.63 


± 





10 


0.62 


0.87 


M 


M (77) 


92 


0.44 


80.668795 


33.491486 


1 


48 


± 





15 


0.59 


± 





12 


0.31 


0.58 


NM 


M (74) 


98 


0.68 


80.722431 


33.509064 


3 


55 


± 





35 


0.50 


± 





09 


0.56 


0.49 


M 


NM (3) 


106 


0.52 


80.751932 


33.496628 


2 


87 


± 





16 


0.53 


± 





06 


0.90 


0.49 


M 


M (74) 


108 


0.46 


80.765402 


33.487503 


3 


13 


± 





34 


0.52 


± 





10 


0.97 


0.56 


? 


NM (45) 


115 


0.44 


80.767404 


33.470600 


2 


90 


± 





24 


0.56 


± 





09 


0.43 


0.38 


M 


- 


123 


0.45 


80.782318 


33.467247 


3 


77 


± 





27 


0.49 


± 





06 


1.56 


0.65 


? 


- 


134 


0.48 


80.744271 


33.400803 


3 


19 


± 





21 


0.57 


± 





08 


0.92 


1.49 


M 


M (79) 


139 


0.40 


80.723816 


33.392075 


2 


77 


± 





19 


0.65 


± 





08 


0.57 


0.10 


M 


M (78) 


140 


0.41 


80.719063 


33.386875 


2 


72 


± 





19 


0.57 


± 





09 


0.75 


0.67 


M 


M (72) 


141 


0.42 


80.717616 


33.384277 


2 


60 


± 





13 


0.57 


± 





06 


0.67 


0.80 


M 


M (78) 


149 


0.77 


80.665351 


33.371723 


2 


81 


± 





07 


0.56 


± 





03 


1.35 


2.43 


M 


M (80) 


168 


0.53 


80.643091 


33.489128 


2 


34 


± 





20 


0.60 


± 





10 


0.66 


0.81 


M 


M (76) 


196 


0.83* 


80.788871 


33.500652 


2 


26 


± 





15 


0.61 


± 





08 


0.73 


1.21 


? 


M (79) 


197 


0.20" 


80.781174 


33.497807 


2 


19 


± 





13 


0.57 


± 





06 


0.96 


0.60 


NM 


NM (0) 


199 


0.73 


80.781643 


33.477089 


2 


82 


± 





24 


0.54 


± 





08 


2.19 


0.98 


? 


- 


209 


0.25* 


80.823050 


33.436195 


2 


43 


± 





30 


0.41 


± 





06 


1.50 


0.46 


NM 


M (80) 


213 


0.53* 


80.798357 


33.402302 


3 


41 


± 





46 


0.51 


± 





11 


0.84 


1.45 


? 


NM (22) 


215 


0.47 


80.789509 


33.406975 


3 


23 


± 





67 


0.41 


± 





10 


0.40 


0.97 


M 


M (81) 


219 


0.48 


80.782304 


33.385307 


3 


81 


± 





36 


0.52 


± 





08 


1.11 


0.78 


M 


M (77) 


222 


0.39 


80.770831 


33.389034 


2 


54 


± 





33 


0.57 


± 





15 


0.21 


1.16 


M 


M (78) 


228 


0.40 


80.741278 


33.368721 


2 


81 


± 





21 


0.62 


± 





09 


1.17 


0.66 


M 


M (77) 


301 


0.39* 


80.831142 


33.452225 


3 


21 


± 





07 


0.55 


± 





03 


1.62 


0.69 


M 


M (78) 


307 


0.42* 


80.854608 


33.435932 


3 


44 


± 





10 


0.57 


± 





04 


0.28 


0.28 


M 


M (79) 


310 


0.20" 


80.847118 


33.419819 


2 


15 


± 





30 


0.45 


± 





07 


0.49 


0.74 


NM 


M (84) 


320 


0.46* 


80.802256 


33.366192 


2 


81 


± 





17 


0.55 


± 





08 


0.25 


0.56 


M 


M (65) 


327 


0.20* 


80.680573 


33.288097 


1 


16 


± 





14 


0.50 


± 





13 


0.68 


0.60 


NM 


M (75) 


329 


0.19* 


80.656141 


33.294849 


1 


15 


± 





20 


0.67 


± 





20 


1.22 


0.83 


NM 


M (76) 


355 


0.20 a 


80.899315 


33.423325 


2 


05 


± 





35 


0.68 


± 





20 


0.09 


0.53 


NM 


NM (48) 


356 


0.20 a 


80.853457 


33.408504 


2 


59 


± 





36 


0.57 


± 





16 


1.18 


0.75 


NM 





t: Cuffey and Shapley (1937) 

i; To estimate E(B — V) values, UBV photometric data has been taken from Sharma et al. (2007) 
*: To estimate E(B — V) values, UBV photometric data has been taken from Masscy ct al. (1995) 
a : The foreground stars with E(B - V)=0.20 mag from ZAMS fitting 
tt; Present Work 
ttt : Dias et al. (2006) 
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Table 10. Mean values of Py and 8y for members and non-members. 



Cluster Id 


P v ± tr(%) 




No of stars 


Members 


NGC 2281 


0.9±0.2 


17±6 


11 


NGC 1664 


1.4±0.2 


172±4 


18 


NGC 1960 


1.2±0.1 


160±3 


12 


Stock 8 


2.5±0.2 


161±4 


13 


NGC 1893 


2.8±0.4 


160±4 


28 


Non-Members 


NGC 2281 


0.9±0.3 


14±8 


3 


NGC 1664 


2.5±0.9 


171±14 


9 


NGC 1960 


1.2±0.4 


155±7 


3 


Stock 8 


1.8±0.3 


161±3 


6 


NGC 1893 


1.6±0.4 


166±4 


11 



